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Abstract Shell convection driven by nuclear burning in a stellar core is a com- 
mon hydrodynamic event in the evolution of many types of stars. We encounter and 
simulate this convection (i) in the helium core of a low-mass red giant during core 
helium flash leading to a dredge-down of protons across an entropy barrier, (ii) in 
a carbon-oxygen core of an intermediate-mass star during core carbon flash, and 
(iii) in the oxygen and carbon burning shell above the silicon-sulfur rich core of a 
massive star prior to supernova explosion. Our results, which were obtained with 
the hydrodynamics code HERAKLES, suggest that both entropy gradients and en- 
tropy barriers are less important for stellar structure than commonly assumed. Our 
simulations further reveal a new dynamic mixing process operating below the base 
of shell convection zones. 



1 Introduction 



Our knowledge of stellar core convection stems from one-dimensional stellar evolu- 
tionary calculations assuming a hydrostatic stellar structure and describing dynamic 
processes (like, e.g., convection) by local and linear theories | Weiss et al. [2004 1. 
This approach is computationally feasible and predicts observables. However, it in- 
volves some degeneracy due to necessary utilization of observationally constrained 
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free parameters |Montalban et al. 2004 1, which may vary during the star's evolution, 
and from star to star. One way to check the ID results is by means of hydrodynamic 
simulations which are based on the solution of the Navier-Stokes equations and are 
essentially parameter free. This approach is computationally demanding, and in our 
case does not yet predict any observables. Nevertheless, our simulations provide 
some important insight into intrinsically multidimensional processes connected to 
core shell convection, like convective overshooting, and mixing across entropy bar- 
riers. In particular, they revealed a potentially very important, new mixing process 
operative below the base of such convection zones. 




Fig. 1 3D volume rendering of the velocity amplitude of shell convection (red/yellow/green) en- 
closed between two stable layers (blue) in the helium core during core helium flash. 



2 Shell convection in stellar cores 



Shell convection in stellar cores occurs in a sandwich-like structure where a dynam- 
ically unstable zone sustained by nuclear burning is enclosed by two stable regions 
above and below (Fig.[T|. Besides mixing of chemical species the convective shell 
redistributes the energy released by nuclear burning and keeps the stellar core in 
quasi-hydrostatic equilibrium. Such convection typically occurs for instance during: 



1. core helium flashes in low-mass red giants | Mocak et al. 2008 2009 2010] 
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2. core carbon flashes of "super- AGB" stars | Siess et al. 2002| 



3. shell nuclear burning in massive stars |Meakin and Arnett 



4. thermal pulses of AGB stars JHerwig 2005 Herwig et al. 



2007] 



2006 



We simulated the first three cases mostly in two but also in three spatial dimen- 
sions (2D and 3D simulations, respectively) using initial core structures of metal- 
rich 1.25M [Mocak et aL||2008) , 9.3 M star piessl[2006| , and 23 M | |Meakin| 
|and A rnett, 2007) stars. None of these simulations led to a violent hydrodynamic 
event. Initially convection occured in regions which are unstable according to the 
Schwarzschild criterion. The convective velocities found in our 3D simulations are 
close to those predicted by mixing-length theory (MLT). However, our simulations 
showed an additional non-radial instability at the boundaries of the convection zone, 
which allowed convection to penetrate inexorably into the adjacent stable layers, 
thereby increasing the width of the convection zone on a dynamical timescale. Con- 
sequently, core convection driven by nuclear burning covers likely larger regions 
than predicted by MLT theory. 



3 Turbulent entrainment and entropy gradients 

Convection is able to generate mixing in neighboring stable layers by convective 
plumes, which can move into these regions due to their momentum. This process 
is typically called overshooting and refers to localized events. However, if the fre- 
quency of these events is high, entropy can change in the affected stable layers. 
Hence, we can speak of penetration |Brummell et al. 2002]. This occurs due to ex- 



tensive heat exchange between the penetrating convective plumes and the stable lay- 
ers. Actually, processes at the edges of convection zones appear to be a combination 
of both overshooting and penetration. We prefer to call them turbulent entrainment 
[Fernando] 1 1991 ) which is well described by the divergence of the buoyancy flux at 
convection boundaries [Meak in~and Arnettl|2007| ]. We find turbulent entrainment to 
operate at convection boundaries in all our hydrodynamic models. We thus conclude 
that convection zones grow in size during dynamic nuclear flashes (nuclear burning 
in semi-degenerate gas) when entropy gradients given by canonical ID stellar calcu- 
lations cannot withstand turbulent entrainment of stable layers. A direct implication 
of this result for the core helium flash is the occurence of a hydrogen injection flash, 
because the entrainment rates found for the convection zone driven by helium burn- 
ing in our models are of the order of meters per second. Thus, the upper boundary of 
the convective shell would reach the overlying hydrogen shell within weeks. This is 
not predicted for solar metallicity stars | Campb ell and Lattanzio| [2008| as our ini- 
tial model. The implications of enlarged convection zones for the core carbon flash 
in intermediate-mass stars remain unexplored, and the case of the oxygen burning 



shell in massive stars was studied by Meakin and Arnett [ 2006 2007 1. 
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4 Hydrogen injection flash and mixing across an entropy barrier 
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Fig. 2 (a) Velocity amplitude | v | (in 1 6 cm s~ 1 ) and (b) hydrogen mass fraction X( 1 H) at the onset 
of the hydrogen injection flash in a meridional plane of a 3D simulation. Panels (c) and (d) show the 
angle-averaged radial distributions of temperature and entropy, respectively, at the same time. The 
two temperature maxima are denoted as l s 'T,„ av and 2"' l T max , and the entropy barrier is marked by 
a vertical arrow. Dashed lines enclose distinguished layers, where "CVZ-1" is the convection zone 
driven by helium burning, "CVZ-2" the convection zone driven by hydrogen burning, and "stable" 
denotes dynamically stable layers. 



Hydrogen injection during the core helium flash is predicted by canonical ID 
stellar evolution calculations only for Pop III and extremely metal-poor stars with 
intrinsic metallicities [Fe/H] < —4]. A similar hydrogen injection phase also occurs: 
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• at the beginning of the thermally pulsing AGB phase of metal-poor intermediate- 
mass stars (^^et^|2^||iiss^^|2^|^^to^^g^ 

• in "Late Hot Flasher" stars experiencing strong mass loss on the RGB | Brown 
[etaLl[200T1[Cassisi et al.| 12003) 



We refer to these events as "dual flashes" | Campbell and Lattanzio 2008 1, since 
they all experience simultaneous hydrogen and helium flashes. These events often 
lead to a "splitting" of the initial convection zone driven by helium burning due to 
dredge-down of material from the above-lying hydrogen-rich envelope and rapid 
CNO burning. 

In order to study this episode by means of hydrodynamic simulations, we de- 
signed a special stellar model from the helium core of a metal-rich 1.25M©, where 
we shifted the hydrogen-rich layers down to the upper edge of the already present 
helium-burning convection zone. This immediately causes hydrogen to be dredged- 
down to the hotter layers of the underlying helium-buming convection zone. The 
results of this first 3D simulation of a dual core flash being driven by hydrogen in- 
jection into hot layers of helium-burning convection zone (Fig.|2]i are the following: 

• dredge-down of protons across an entropy barrier between the helium-rich layers 
and the hydrogen shell 

• subtle retreat of the helium-burning convection zone (CVZ-1 in Fig. [2} to smaller 
radii, and the appearance of a hydrogen-burning convection zone (CVZ-2 in 
Fig. |2]i just above it 

• non-existence of stable radiative layer preventing mixing between the two con- 
vection zones, which is in contradiction to ID stellar evolutionary calculations 

• appearance of the 2 nc ^ temperature peak at the base of the hydrogen-burning 
convection zone, in agreement with ID stellar evolutionary calculations 

We think these results should be qualitatively similar to all dual flash events and to 
thermal pulses occurring during AGB phase of stars. Hydrodynamic simulations of 
a thermal pulse with hydrogen mixing into the helium-burning convection zone has 
already led to improvements in our understanding of Sakurai's object [Herwig et al.| 
|20TT) . 



5 A new dynamic mixing process below the base of shell 
convection zones 

In our 2D and 3D hydrodynamic models of the core helium and carbon flash, we 
discovered unreported dynamic mixing process operating at the base of convection 
zones. The mixing manifests itself by cold and over-dense blobs sinking in the di- 
rection of gravity, leaving traces of material with higher mean molecular weight 
fl (Fig.|3j. Blobs originate from a layer situated just below the lower edge of the 
convection zone. We do not see any mixing in our oxygen burning shell models. 
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Fig. 3 Relative fluctuations of the mean molecular weight defined as A[l/ [l = (fl — (fi)e) I 
(taken from 2D simulations) below the base of shell convection zones during the (a) core helium 
flash, (b) core carbon flash, and (c) oxygen burning shell shown together with the radial distribu- 
tions of the mean molecular weight fl (solid) and the square of the Brunt- Vaisala frequency N 2 . 
The horizontal dashed line corresponds to N 2 = 0, and ()g denotes the horizontal average at a given 
radius. 

The driving mechanism of the mixing remains unclear. Possible explanations for 
the appearance of the sinking cold and dense blobs range from: 

• the presence of a steep negative mean molecular gradient < with = 
dlnjj,/dlnP (i.e., the molecular weight ji decreases in direction of gravity) desta- 
bilizing the layers where the mixing starts 

• a strong shear creating peculiar turbulence at the convection boundary due to 
trapped gravity waves excited by convection in the zone above and indicated by 
a large value of the square of the Brunt-Vaisala frequency N 2 

• a weak dynamic stability of the layers below the convection zone indicated by 
a relatively small positive value of N 2 , which is larger in case of the oxygen 
burning shell model, and hence mixing is effectively suppressed 

• a "turbulence" imitating heat transfer or grid scale diffusion removing heat from 
patches of gas at the convection zone boundary, hence making some blobs colder 
and heavy 
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6 Summary 

We investigated convection driven by shell nuclear burning during core helium flash 
in a low-mass star, core carbon flash in an intermediate-mass star, and in the oxygen 
burning shell in a massive star using data from 2D and 3D hydrodynamic simula- 
tions. We find that all convection zones are growing on a dynamic timescale due 
to turbulent entrainment. This either implies that the convection zones become en- 
larged or that the physical conditions at the convection zone boundaries inherited 
from the initial ID stellar evolutionary models are inconsistent with real conditions 
in stars. The enlargement of the helium-burning shell convection zone during the 
core helium flash may lead to a dual core flash by injection of hydrogen into deeper 
hot layers of the helium core. With our specifically designed initial model for a hy- 
drogen injection flash we show that such an event is possible and the entropy barrier 
does not prevent mixing of material from the hydrogen shell into the underlying hot 
helium-burning convection zone. These findings imply that the stabilizing entropy 
gradients are less important in stars during flash events. In the hydrogen injection 
flash simulation, we observe the appearance of a second hydrogen-burning shell re- 
siding above the already existing helium-burning convection zone. Radiative stable 
layers separating the two convection shells are not observed. We think, that these 
finding should be generally valid for all dual flashes. We also found a new dynamic 
mixing process below the base of shell convection zones manifesting itself by cold 
and dense blobs sinking in direction of gravity creating eventually finger-like struc- 
tures. The origin of this mixing is still under investigation. 
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